Conductance of Quantum Impurity Models from Quantum Monte Carlo 
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The conductance of two Anderson impurity models, one with two-fold and another with four- fold degeneracy, 
representing two types of quantum dots, is calculated using a world-line quantum Monte Carlo (QMC) method. 
Extrapolation of the imaginary time QMC data to zero frequency yields the linear conductance, which is then 
compared to numerical renormalization group results in order to assess its accuracy. We find that the method 
gives excellent results at low temperature (T < Tk) throughout the mixed valence and Kondo regimes, but it is 
unreliable for higher temperature. 
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Quantum dots provide a highly controlled and tunable 
way to study a range of quantum many-body physics: var- 
ious quantum impurity models and their associated Kondo 
effects,'"^ tunneling with dissipation,^ and Luttinger liquid 
effects, ^-^ to name a few. The crucial experimental observ- 
able in these situations is the conductance; thus, calculat- 
ing the conductance is a key task for both analytic and nu- 
merical approaches. Numerical methods have indeed been 
developed, '^'^ with remarkable agreement for small systems 
between theory and experiment.'"* But these methods scale 
poorly for the larger, more complex multi-dot systems'^''^ that 
are currently of great interest. Here we implement and test a 
way to calculate the conductance from a path-integral quan- 
tum Monte Carlo (QMC) calculation. While it yields less in- 
formation than NRG in simple systems (e.g. a single quantum 
dot), the method should scale readily to more complicated sys- 
tems. Results for two Anderson-type impurity models show 
that the method works very well at low temperature. 

For calculations of the conductance in simple quantum dot 
systems, the most accurate results are obtained using the nu- 
merical renormalization group (NRG) method. '^'^^ NRG be- 
comes slow and even impractical, however, if there are many 
leads, a many-fold degeneracy, or more than a few interact- 
ing sites. In such situations, the world-line quantum Monte 
Carlo (QMC) method could be a valuable alternative since it 
scales nicely as the problem size increases. However, QMC is 
formulated in imaginary time rather than real time: to extract 
dynamic properties one must transform from imaginary back 
to real time. The statistical error in the QMC data makes this 
an ill-posed problem, for which various extrapolation and con- 
tinuation methods have been developed.'^ To obtain the con- 
ductance of interest here, we extrapolate to zero frequency the 
appropriate correlation function evaluated using QMC at the 
imaginary Matsubara frequencies^^ This has been used, for 
instance, to study a one-dimensional Hubbard chain coupled 
to non-interacting leads in the absence of the Kondo effect. ^^ 

The aim of this paper is to test the validity of the extrapola- 
tion method for Anderson impurity models in both the mixed 
valence and Kondo regimes. We study the linear conductance 
using QMC in two models: a single impurity Anderson model 
with either two-fold or four-fold degeneracy. The standard 
two-fold degenerate model is a simplified representation of a 
single GaAs quantum dot connected to leads.' The four- fold 
degenerate model represents a quantum dot in a carbon nan- 



FIG. 1: (color online) The ID infinite tight-binding chain, where the 
O"* site is the impurity site (quantum dot). 



otube in which there is an additional orbital degeneracy from 
the helicity of the states.'^ This orbital degeneracy is present 
in both the discrete states in the dot and the extended states in 
the carbon nanotube leads. 

Consider a model, then, in which a single level with 
Coulomb repulsion U represents the quantum dot (which we 
also refer to as the impurity site) and is coupled to two non- 
interacting bands, left (L) and right {R). The degeneracy of 
both the discrete level and the free electrons is M; we will 
consider the two cases M — 2 (standard single-level Anderson 
model) and M ~ i (both spin and orbital degeneracy). The 
Hamiltonian is 
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where the electron number operator for the impurity site is 
N — J2a=i ^\^<y The energy in the bands is such that —D < 
Ck < D where D is the half bandwidth, and we assume a 
flat density of states, p = 1/2D. The hybridization of the 
impurity to each lead is given by V^ which yields a level width 
Ta='^L,a + TR,a With T L,cr ^'^ R,a = T^V^ p. In terms of the 
gate voltage Ng, the energy level of the dot is explicitly given 
by Cd = U{l — 2Ng)/2. Finally, in the absence of any orbital 
degeneracy, the degeneracy of the d level is simply given by 
spin, a =t or |. 

Method — A new basis for the two non-interacting bands can 
be independently constructed by starting from the localized 
impurity state. In this way the model is mapped to a one- 
dimensional (ID) infinite tight-binding chain,'" as shown in 
Fig.[T] We use a large chain (^ 10^ sites) in order that its 
finite size is irrelevant for the physics of interest. Then, in 
order to make the computation time manageable, logarithmic 
blocking of the energy levels'" is used to reduce the number of 



effective sites. In this work, the logarithmic blocking factor is 
A = 2.5 (the number of effective sites is ^ 61). We use a form 
of blocking^^ which avoids A-dependent corrections^! to the 
low-energy scales [i.e. Tk(A)]. We solve the resulting prob- 
lem using the world-line quantum Monte Carlo method with 
a directed-loop cluster algorithm.i^^ The Trotter number N 
is choosen such that e = f3/N^0.1/D. 

To find the conductance, we proceed following the method 
of Syljuasen in Ref.|22 which is itself closely related to sev- 
eral other approaches. '^^' The conductance at the (imagi- 
nary) Matsubara frequencies, g{iujn) with uj„ — 2TrnT, is re- 
lated in linear response to the current-current correlation func- 
tion in the usual way. For a one-dimensional system with open 
boundary conditions, current continuity can be used''*-^**-^' to 
express .g(ia;„) in terms of charge correlations (polarizability). 



7(ia;„) = — 



drcos(w„r)(P,(T)Py(0)) (2) 



where Py is the sum of the electron charge density operators 
to the right of y, Py = X]t,'>w "-t,' ■ Thus the time derivative 
of {Py) is the current through the bond between sites y — 1 
and y. We calculate g{iuin) for n > from the world-line 
QMC data in this way. Not all combinations of x and y can 
be used in Eq. ^ because the system is not a physical chain 
but only effectively mapped to a chain. Notice that the current 
through the four bonds closest to the impurity site (labeled 0) 
correspond to the physical current. Therefore, x and y must 
be chosen from among { — 1,0, 1,2}. In addition, left-right 
symmetry reduces the number of independent combinations. 
In our calculation, we choose three cases for x and y: (0, 1), 
(0,0), and (-1,0). 

The hnear conductance G is obtained by extrapolating 
g{iujn) to zero frequency, G — lim^^^o g{iujn)- We carry 
out this extrapolation as follows. First, we try to fit the data 
at the four or five lowest Matsubara frequencies [g(ia;„) for 
n = 1, ..., 4 or 1, ..., 5] to a linear or quadratic polynomial. If 
this method yields a good fit, we simply extrapolate the data 
by using the polynomial. If neither polynomial fit is good, 
the data at the first 14 lowest Matsubara frequencies are fit by 
using a series of rational polynomial functions of different de- 
gree [p/q] (e.g. p for the numerator, q for the denominator, 
p = 1 for a constant, p = 2 for linear function, etc.) as de- 
scribed in Ref.i22| We use allp and q such that 5 < p+q < 10 
and p,q > 2 but exclude cases in which spurious poles ap- 
pears. The final extrapolated value is the average of the results 
for these different forms, and the error bar at zero frequency is 
the maximum spread, which is larger than the error bar of any 
single [p/q] extrapolation. To justify this method, we check 
that three conditions are met. (i) The data for all the combi- 
nations of X and y must extrapolate to nearly the same value 
(the current through different bonds at non-zero frequency can 
be different, but current continuity requires that at zero fre- 
quency the current through all bonds be the same), (ii) The 
data should fit well to most of the functional forms of degree 
[p/q] (we cannot exclude too many cases), (iii) Finally, the 
conductance should have a small error bar (a large error bar 
shows that the extrapolation is model dependent). 
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FIG. 2: (color online) Conductance through a single-level Anderson 
model without orbital degeneracy as a function of gate voltage: QMC 
result (symbols) compared to NRG calculation'^ (lines). Data for 
four temperatures are shown: T ~ 0.11 Tk (brown, dot-dashed, /3 = 
98.3), T ~ 0.24 Tk (blue, upper solid, ^ = 43.7), T ~ 1.2 Tk (red, 
dashed, 13 = 8.6), and T ~ 6.1 Tk (black, lower solid, /3 = 1.7). For 
T ~ 6.1 Tk, the black stars are for {x, y) — (0, 1) while the black 
circles are for (0,0). Tk denotes the Kondo temperature found by 
NRC''' at the particle-hole symmetric point (— 6^/(7 = 0.5). Note 
the high accuracy of the QMC result as long as T<Tk. 



Conductance without orbital degeneracy — We first con- 
sider the standard single-level Anderson model, M = 2 in 
Eq. ([T]). We compare the conductance obtained by our QMC 
calculation to that from the numerical renormalization group 
(NRG) calculation of Ref.Uj Tsee their Fig. 2(a)]. The param- 
eters are D = 100, T = 1.0, and U = Svr. The NRG valuell 
for the Kondo temperature at the particle-hole symmetry point 
(— e/C/ = 0.5), which we denote Tk throughout, is Tk — 0.1. 

Fig.|2] compares our calculation of the conductance as a 
function of gate voltage to the NRG results '"* for several tem- 
peratures. The QMC results are in excellent agreement with 
the NRG results for T < TKfor all values of the gate voltage — 
that is, in both the mixed valance and Kondo regimes. For T 
slightly larger than Tk, agreement is good; in contrast, note 
that there is a substantial error in the extrapolated conductance 
value for larger T. 

Some examples of the extrapolations used to obtain the con- 
ductance shown in Fig. |2] are given in Figs.[3]|5j moving from 
lower to higher temperature. Fig. [3] shows four examples of 
the conductance at imaginary frequency, g(za;„), for T < Tk- 
Examples of a linear fit [panel (a)], a quadratic fit [panel (c)], 
and rational polynomial fits [panels (b) and (d)] are shown. In 
the mixed valance regime, —ed/U < 0.1 or > 0.9, a linear or 
quadratic polynomial works well, and the three curves for dif- 
ferent (x, y) all extrapolate to nearly the same value, leading 
to a small error bar In the Kondo regime, 0. 1 < — e^/C/ < 0.9, 
the linear or quadratic polynomial does not fit well, but the 
QMC data can be fit to a series of rational polynomials as dis- 
cussed above. Almost all values of [p/q\ work well, and the 
three sets of g{iLdn) for different {x, y) extrapolate to nearly 
the same value, leading to a small error bar In this tempera- 
ture regime, then, the extrapolation is straight forward and the 
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FIG. 3: (color online) Conductance at Matsubara frequencies at low 
temperature (symbols) for the single-level Anderson model without 
orbital degeneracy and the corresponding fits used to extrapolate to 
zero frequency (lines). The values of T and — e/C/ are (a) 0.11 Tk, 
-0.1; (b) 0.11 Tk, 0.5; (c) 0.24 Tk, 0.1; (d) 0.24 Tk, 0.3. Points 
for three choices of {x, y) are shown: (0, 1) red triangles, (0, 0) blue 
squares, and (—1,0) green circles. A good quality extrapolation is 
obtained in all cases. 



FIG. 5: (color online) Conductance at Matsubara frequencies for 
high temperature, T ~ G.ITk, in the absence of orbital degeneracy 
(symbols) and the corresponding fits used to extrapolate to zero fre- 
quency (lines), (a) — e/(7 = — 0.1; (b) — e/f/ = 0.5. Points for three 
choices of {x, y) are shown: (0, 1) red triangles, (0, 0) blue squares, 
and ( — 1, 0) green circles. The black stars are the NRG data. The 
accuracy of the extrapolation is poor in all of these cases. 
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FIG. 4: (color online) Conductance at Matsubara frequencies for 
T ~ 1.2 Tk in the absence of orbital degeneracy (symbols) and 
the corresponding fits used to extrapolate to zero frequency (lines). 
The values of — e/t/ are (a) 0.1 (mixed- valence) and (b) 0.5 (Kondo 
regime). Points for three choices of [x, y) are shown; (0, 1) red tri- 
angles, (0, 0) blue squares, and ( — 1,0) green circles. Extrapolation 
using {x, y) = (0, 1) is accurate. 



agreement with the NRG result is excellent. 

For T ^ Tk, two examples of the conductance function 
g{iujn) are shown in Fig.|4] For the mixed valance regime 
[Fig.|4ja)], a quadratic polynomial works well for {x,y) — 
(0,0) and (—1,0), and rational polynomials are used for 
(0,1). All three combinations extrapolate to nearly the same 
value, so the result is accurate. In the Kondo regime [panel 
(b),— Erf/C/ = 0.5], g{iuJn) for (0, 1) can be fit with rational 
polynomials. However, for both other cases, {x, y) = (0, 0) 
and (—1,0), there is a small wiggle near cu — 2ttTk in the 
imaginary frequency conductance function g{iujn), showing 



that there is important structure below that frequency. Since 
there is only one data point below uj = 2ttTk, the extrapola- 
tion is unreliable. Thus, we do not use the data when struc- 
ture appears at a frequency below which there are only a few 
data points. The conductance in the Kondo regime for this 
temperature is based only on (a;, y) = (0, 1); nonetheless, the 
agreement with the NRG result is good. 

Finally, for T > Tk (Fig.©, the functions g{iuJn) for the 
three combinations of {x, y) do not extrapolate to the same 
zero-frequency value. Notice also that the conductance ob- 
tained in the mixed-valence regime (the gate voltage at which 
the conductance peaks for this temperature) has a large error 
bar. For the cases (x, y) = (0, 0) and (-1, 0), the QMC data 
can be fit with a rational polynomial, but the extrapolated re- 
sult disagrees substantially with NRG. For the case (0, 1), the 
average value of G from QMC roughly follows the NRG re- 
sult (Fig.|2]), but the large error bar in most cases indicates that 
the result has little meaning. Thus, the QMC extrapolation 
method is unreliable for T substantially larger than Tk- 

Conductance with orbital degeneracy — We now turn to 
considering an Anderson model in which all the states, both 
those in the dot and in the leads, have an orbital degeneracy 
in addition to spin degeneracy: i\f = 4 in Eq. ([T]). This sit- 
uation arises, for instance, in carbon nanotube quantum dots 
connected to carbon nanotube Ieadsi2rdi2i To assess the qual- 
ity of our QMC results, we compare with the NRG results of 
Ref. 12 (see their Fig. 16). The parameters we use are D — 30, 
U = O.ID = 3, ri,2 = 0.003 ttA and r3,4 = 0.002 ttD. 
At the particle-hole symmetric point where the Kondo tem- 
perature is a minimum, the NRG estimation'^ for Tk yields 
Tk ~ 0.0042.. 

Fig.|6] compares our calculation of the conductance as a 
function of gate voltage to the NRG'^ results. For T < Tk, the 
QMC and NRG results are in very good agreement through- 
out both the mixed valence and Kondo regimes. For T >Tk 



4 

3.5 

3 

_2.5 

& 

(3 
1.5 

1 

0.5 





-•-NRG (T=0.30Tk) 




■. QMC (T=0.30Tk) 


ft 


-»-NRG (T=0.93Tk) 


h*' \ 


n QMC (T=0.93 Tk) 


-0- NRG (T=2.7 Tk) 


« 


tV QMC easel {T=2.7 Tk) 


r-A -- 


QMC case2(T=2.7 Tk) 



-0.05 




0.35 0.4 



FIG. 6: (color online) Conductance in four-fold degenerate model as 
a function of gate voltage: QMC results (symbols) compared with 
NRG calculations'' (lines with symbols). Results for three tempera- 
tures are shown: T~ 0.30 Tk (blue circles and solid line, /3 = 79.4), 
0.93 Tk (red squares and dashed line, /3 = 25.6), and 2.7 Tk (green 
triangles or stars and dotted line, /3 — 8.77). For the highest tem- 
perature (T ~ 2.7 Tk), the QMC data labeled case 1 are based on 
{x, y) — (0, 0) and (—1, 0) while those for case 2 use (0, 1). Tk here 
denotes the Kondo temperature found by NRG*" at the particle-hole 
symmetric point. The good agreement of the QMC data with the 
NRG results illustrates the value of the QMC approach, though note 
the growing error bar when T > Tk- 



(T ~ 2.7 Tk), the QMC conductance roughly follows the 
NRG result but does not accurately agree with it. In addition, 
a large error bar is encountered at the highest temperature, 
showing that, as in the doubly degenerate case, the extrapola- 
tion is not reliable for these temperatures. 

Four examples of the extrapolation from the imaginary fre- 
quency conductance function, g{iuj„), are shown in Fig.|7] At 
low temperature, panel (a), the extrapolation is good and con- 
sistent for all three values of {x, y) using the rational polyno- 
mial fit. Near the particle-hole symmetry point and for T ~ Tk 
[panel (b)], the case with (a;, y) — (0, 1) fits nicely to a ratio- 
nal polynomial and the extrapolated value agrees with NRG. 
For the other two curves ((x, y) = (0, 0) and (—1, 0)), a small 
wiggle appears near uj ~ 27rTK, as in the case without orbital 
degeneracy (Fig.|4]i, making extrapolation difficult. For larger 
T, panels (c) and (d), although the QMC data for two cases 
((x, y) = (0, 0) and (—1, 0)) can be fit to rational polynomi- 



als and yield an estimated conductance with small error bar, 
the value does not agree accurately with the NRG result. The 
(x, y) = (0, 1) yields a large estimated error. Therefore, as we 
saw in the case without orbital degeneracy, when the temper- 
ature become large, the QMC method becomes inaccurate. 

In summary, we developed and tested a method to ob- 
tain the linear conductance by extrapolating from QMC data. 
By studying two cases for which NRG results exist in the 
literature r^iil we demonstrated the accuracy of the extrapo- 
lation technique as long as the temperature is not too high, 
T < Tk (where Tk denotes the Kondo temperature at the 
particle-hole symmetric point). We expect that this technique 
will be useful for finding the conductance of more complex 
quantum dot and/or impurity systems, such as three and four 
quantum dot structuresJ^ii^ 

We thank R Bokes and J. Shumway for helpful discus- 
sions. This work was supported in part by the U.S. NSF Grant 
No.DMR-0506953. 
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FIG. 7: (color online) Imaginary frequency conductance function for 
single impurity Anderson model with orbital degeneracy {M — 4). 
The values of T and -e/U are (a) 0.30 Tk, 0.06, (b) 0.93 Tk, 0.15 
(near particle-hole symmetry), (c) 2.7rK, 0.15, (d) 2.7rK, 0.02. 
Points for three choices of (x, y) are shown: (0, 1) red triangles, 
(0,0) blue squares, and (—1,0) green circles. The black stars are 
the NRG data. The extrapolation is successful at low temperature but 
becomes increasingly problematic at higher temperature, T>Tk. 
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